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1.  INTRODUCTION 


A  problem  of  general  interest  is  to  be  able  to  use  a  set  of  L  beams,  steered  and  focused 
respectively  at  the  bearing-range  pairs  [62,  R2],  — ,  to  approximate  a  beam 

steered  and  focused  at  a  given  angle-range  pair  A  problem  of  particular  interest  is 

to  be  able  to  use  L  beams  focused  at  infinite  range  (plane  wave)  to  approximate  the  beam 
steered  and  focused  at  the  angle-range  pair  [^Oj  -Ro]- 

To  do  this,  a  set  of  L  appropriate  coefficients,  which  are  complex  in  general,  must  be 
computed  and  then  used  to  scale  and  combine  the  beams  to  generate  an  approximation  of  the 
beam  in  question.  These  coefficients  are  derived  in  this  report  via  mean  square  error  (MSE) 
minimization.  Also,  the  number  of  coefficients  needed  for  a  given  error  value  is  examined  as 
a  function  of  frequency  and  focusing  range.  The  beam  interpolation  scheme  is  applied  to  the 
K-u>  beeimformer. 

Throughout,  vectors  are  denoted  by  bold,  lowercase  letters  and  matrices  by  bold,  upper¬ 
case  bold  letters.  All  others  are  scalars.  The  complex  notation  j  (where  j  —  \f^)  is  used. 
The  superscript  ('^)  is  used  to  denote  vector  and  matrix  transpose. 
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2.  GENERAL  OBSERVATIONS 


Almost  certainly,  the  steering  angle  9q  of  the  approximated  beam  must  lie  inside  the 
bearing  span  of  approximating  beams  [di,  9l]. 

The  choice  of  beam  bearing  spacing  for  the  set  of  reference  beams  is  critical  for  the 
proper  formation  of  the  desired  beam  steered/focused  at  a  given  bearing/range.  If  the  ref¬ 
erence  beams  are  spaced  too  far  apart  in  bearing,  then  information  is  lost  and  the  desired 
beam  will  not  be  properly  formed.  If,  on  the  other  hand,  the  reference  beams  are  spaced 
too  closely  together  in  bearing,  then  the  information  provided  by  each  beam  becomes  highly 
correlated.  This  usually  results  in  large  positive  and  negative  complex  coefficients,  leading 
to  numerical  instability. 

The  ideal  spacing  for  a  given  frequency  occurs  when  the  beam  responses  are  orthogonal. 
For  a  fixed  beam  spacing  this  condition  can  be  satisfied  for  only  one  frequency.  For  the 
other  frequencies,  the  reference  beam  set  will  either  be  under-sampled  or  over-sampled.  For 
the  case  of  a  K-u  beamformer,  the  beams  in  K  (bearing)  space  are  always  orthogonal  for 
all  frequencies,  which  follows  from  the  Fourier  transform  property.  This  makes  the  K-uj 
beamformer  an  ideal  choice  for  implementing  this  technique.  This  is  especially  true  when 
the  reference  beam  set  is  formed  at  infinite  range. 
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3.  COEFFICffiNT  DERIVATION 


3.1  FORMULATION  OF  THE  MEAN  SQUARE  ERROR  (MSE) 

A  steered  and  focused  narrowband  beampattern  for  an  iV-channel  line  array  is  given  by 

N 

1=1 

The  time  delay  Tj  {6,  R)  is  the  time  required  for  a  spherical  wavefront  of  radius  R  to  travel 
at  a  speed  c  from  the  center  of  the  array  (at  coordinate  0)  to  element  i  (at  coordinate  Xi). 
Denoting  the  distance  between  the  wavefront  source  to  element  i  by  Ri  (see  figure  1),  one 
has 

ni0,R)  =  ^^. 

c 


Figure  1 .  Time  Delay  Computation  for  a  Curved  Wavefront 


By  using  the  law  of  cosines,  Ri  can  be  expressed  as  a  function  of  Xi  and  the  source  range 
R.  Completing  the  square  and  factoring  (i?  —  Xi  cos  9)  out,  the  time  delay  becomes 


Ti{e,R)  = 


_  1 


R-iR-XiCOs9)Jl+{^^) 


xf  sin^  6 


c 

Xi  cos  6 

c  2c(R  —  Xi  cos  0)  ’ 

curvature 


R'^Xi 


(2) 


where 
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ri{e,R) 

ri{Oo,Ro) 

Xi 

c 

f 

Wi 


Target  time  delay  for  a  target  at  bearing  6  and  range  R, 

Beamformer  time  delay  for  a  steering  angle  and  focus  range  Rq, 
Coordinate  of  element  i  with  respect  to  a  given  reference, 

Speed  of  sound  (beamformer  and  actual  speeds  of  sound  are  assumed  equal). 
Target  frequency,  and 

Shading  weight  applied  to  channel  i,  normalized  so  that  =  1- 


The  approximation  in  equation  (2)  is  the  basic  yT+a  »  1  +  f  for  a  «  1,  which  is  true  in 
this  case  for  ranges  large  enough  when  compared  to  array  dimensions. 


The  generally  complex  coefficients  that  will  be  used  to  combine  the  T-beams  to  approxi¬ 
mate  a  single  beam,  can  be  derived  by  minimizing  the  MSE  between  the  actual  beam  desired, 
steered  to  and  focused  at  Rq,  and  the  weighted  sum  of  L  beams,  steered  to  [^i,  ^2,  ■■■  -,^1] 
and  focused  at  [i?i,  i^2^  •  •  • , -Rl],  respectively.  However,  to  simplify  derivation  and  imple¬ 
mentation,  minimization  of  the  sum  over  elements  of  the  individual  MSEs  for  the  respective 
actual  versus  interpolated  beams  will  instead  be  performed.  Denoting  the  interpolation  coef¬ 
ficients  by  [ai  +jPi,  0:2  +jP2)  ••• ,  cxl  +  JRl],  the  sum  over  elements  of  the  individual  MSEs, 
which  must  be  minimized  over  the  coefficient  ap  j0p,  can  be  expressed  by 


JV 

^  =  E 


i=l 


^.^2wflTiie,R)-Ti{eo,Ro)]  _  + 

P=1 


(3) 


It  is  assumed  in  equation  (3)  that  the  sensor  elements  of  the  L  beams  —which  are  weighted 
and  summed  to  produce  the  new  beam—  are  unshaded  and,  hence,  the  resultant  coefBcients 
will  include  (i.e.,  perform),  the  shading,  in  addition  to  the  focusing-steering  operations.  For 
the  case  of  obtaining  coefficients  that  do  not  include  shading,  a  slighly  different  vector/matrix 
definition  must  be  used  in  the  solution  equation  (9),  and  this  definition  is  given  in  appendix  A. 


Factoring  out  the  inner  terms  of  equation  (3)  and  replacing 

the  magnitude  operator  by  the  equivalent  product  of  the  inner  quantity  and  its  conjugate, 
and  using  the  substitution  Vi{p)  =  27r/[ri(<9o,  Rq)  -  R^)],  one  obtains 


E 


N 

E 

t=l 

N 


E 


L 

L 

Wi  -  E(% 

P=1 

q=l 

L  L 

-  Wi  J2iap  + 

9=1  p=l 


L  L 

+EE(‘^.  +  3Pp){.ocq  - 

p=lq=l 


(4) 
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Define  the  L  x  1  vectors  a,  13,  c,  and  s  as  follows: 


a  =  ,aL]^ , 

P  =  ,01^, 

■  N  N  N 

c  =  '^WiCOSVi{l),y^WiCOsvi{2),---  ,y^Wi  cos  Vi{L) 


s  =  Wj  sint;i(l),  Wj  sint;i(2),  sin  Vi{L) 

A=1  i=l  i=l 

In  addition,  define  the  Lx  L  matrices  C  and  S  so  that 


[C]p5  =  Scos[^;i(g) -'yi(p)]  =  [Cj^p, 

(6) 

[S]p9  =  I^sin[r;i(9)-i;i(p)]  =  -[8],^. 

Z=1 

Also,  for  general  L  x  1  vectors  a,  b  with  p**'  component  notations  [a]p  and  [b]p,  respectively, 
and  the  Lx  L  matrix  A  with  p?***  component  notation  Ang,  one  has 


E 

P=1 


=  a^b, 


EENp[bI,[A)„  =  a-^Ab. 

p=l9=l 

Finally,  using  the  Euler  identity  =  COS7  +  jsin7  to  convert  all  complex  exponentials  to 

cosines  and  sines,  the  scalar  E  in  equation  (4)  can  be  expressed  in  vector  notation  as 

E  =  ^w?-  2a'^c  +  2/3'^s  +  a^Ca  +  p^Cp  -  a^Sp  +  p'^Sa.  (7) 

i=l 

3.2  MINIMIZATION  OF  MSE  AND  DERIVATION  OF  COEFFICIENTS 

To  minimize  F  as  a  function  of  a  and  /3,  the  derivative  of  E  with  respect  to  a  and  P 
is  computed  and  set  to  zero.  The  vector  and  matrix  gradients  identities  of  Scharf  [1],  which 
will  be  used  in  minimization,  are  repeated  here  for  convenience: 

?i:b'^x  =  Jrx'^b  =  b. 


l^x^Ax  =  2Ax, 

where  A,  an  L  x  L  matrix,  and  b,  an  L  x  1  vector,  are  independent  of  x,  an  L  x  1  vector. 
Setting  -^E  =  0  and  ^E  =  0  gives 


4.  APPLICATION  TO  LINE  ARRAY 


4.1  K-lo  background 

The  digital  K-uj^  or  wavenumber- frequency,  field  is  generated  by  way  of  a  two-dimensional 
Fast  Fourier  Transform  (FFT)  over  time  and  space.  The  temporal  FFT  is  applied  on  an  Nf^- 
block  of  the  sampled  channel  output,  where  A^fft  is  a  radix-two  number.  Usually,  the  method 
of  overlap-and-save  is  used  for  the  time-series  data.  The  spatial  FFT,  however,  is  applied 
along  elements,  where  zero-padding  is  usually  used  to  increase  the  number  of  elements  N  to 
the  radix- two  number  Mfft  for  FFT  application.  Care  must  be  taken,  however,  when  com¬ 
puting  the  matrices  C  and  S  in  equation  (6)  after  zero-padding  the  array  elements,  to  set 
the  matrix  entries  corresponding  to  the  padded  zeros  to  zero. 

The  relationship  between  K-uj  beamforming,  which  is  the  two-dimensionald  discrete 
fourier  transform  (DFT)  over  space  (elements)  and  time  (element  outputs),  and  time-delay 
beamforming  arises  from  the  similarity  between  the  types  of  beamforming  and  the  DFT 
equations.  For  a  line  array  with  N  elements  indexed  equi-spaced  at  distance  d.  and  shaded 
by  Wi^  the  planewave  beamformer  output  for  an  incoming  source  signal  Xi{t)  measured  at 
element  i  is  expressed  by 

N 

'Y^wiXi[t  -  Ti{e)],  (10) 

i=l 

where  the  beamformer  time-delay  for  a  planewave  steered  to  angle  0  is 

Ti{6)  ~~(^^^~^)^  i  =  (11) 

- V - ^ 

Xi 

Applying,  first,  the  temporal  Fourier  transform  to  the  beamformer  in  equation  (10), 
substituting  for  Ti{9)  in  equation  (11),  and  shifting  the  index  i  from  the  range  [1,A'’]  to 
[0,  iV  —  1]  yields 

g-j27r/^f  cose  ^  (12) 

i=0 

where  /  is  the  frequency  variable  and  Xi{f)  is  the  temporal  Fourier  transform  of  the  signal 

Xi{t). 


To  implement  the  sum  over  elements  in  equation  (12)  using  the  spatial  Mf^-length  FFT, 
with  Mfft  —  N  zero  padding,  the  following  change  of  variables  must  be  performed: 

f-cosO  =  -—,=^ 

C  Mfft 


,  _  MstdcosB  j. 

^  ~~  I 
c 
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which,  for  a  target  at  the  array  design  frequency  (A  =  2d  =  c/f),  or  any  other  smaller  fre¬ 
quency,  will  lead  to  &  k  e  [-^,  ^].  Negative  values  for  k  indicate  0  >  90°,  while  positive 
ones  indicate  6>  <  90°;  A:  =  0  indicates  a  broadside  target.  Since,  however,  the  FFT  considers 
the  integer  range  [0,  - 1],  the  9  >  90°  values  will  appear  in  the  range  Mfft],  while 

the  9  <  90°  values  will  appear  in  the  range  [0,Mfft/2  -  1]. 


Going  back  to  the  FFT  implementation  of  K-u,  after  substituting  N  with  the  radix-2 
Mfft,  WiXiif)  with  y(z)  for  convenience,  and  the  values  of  A:  from  equation  (13),  equation  12 
representing  the  K-u  field  X{k^  /)  becomes 

(14) 

i=0 

To  implement  equation  (14)  via  a  forward  FFT,  the  index  i  must  first  be  flipped  from 
[0,Mfft  -  1]  to  [-(Mfft  -  1),0].  Then,  performing  the  index  substitution  i  i  -  (Mfft  -  1) 
and  factoring  the  non-index-dependent  phase  out  of  the  summation,  equation  (14)  becomes 

X{k,f)=e^^'^at  ^  y(^Mfit  -  1  -  ,  (15) 

i=:0 

which  may  be  implemented  via  a  forward  FFT.  The  term  y{Ms^  —  1  —  i)  represents  a  flipped 
version  of  y{i),  which  is  easily  obtained  by  applying  the  spatial  FFT  with  the  element  order 
reversed. 


For  the  frequency  axis,  the  usual  bin-to-frequency  formula  is  used: 


where  /  is  the  frequency  of  interest  corresponding  to  bin  i,  Fg  is  the  temporal  sampling 
frequency,  fVfft  is  the  FFT  length,  and  Z  =  0, 1,  •  •  • ,  iVgt  —  1  is  the  frequency  bin  index. 

Each  point  of  the  K-u  field  defined  by  its  coordinates,  frequency  bin  =  I  and  wavenumber 
bin  —  k,  represents  a  frequency  filter  that  allows  targets  at  frequency  corresponding  to  bin 
i  and  a  spatial  filter  that  allows  targets  at  bearing  obtained  by  substituting  k  in  equation  (13). 

The  location  of  the  spatial  filter  is  a  function  of  the  frequency  of  interest.  For  different 
frequencies  /  =  c/A,  a  given  maximum  response  axis  (MRA)  corresponds  to  different  k  values 
defined  by  equation  (13).  More  specifically,  equation  (13)  prescribes  a  linear  relationship 
between  the  frequency  (c/ A)  and  wavenumber  k.  Thus,  a  broadband  target  with  bearing  9 
appears  in  the  K-uj  spectrum  on  a  straight  line  defined  by 
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(16) 


fc  = 


Mfft^COSg  y*  0^  <  0  <  90° 

Msihsilf  +  Mfft  90°  <e<  180° 


(A:  =  0,L---,Mfft/2-l), 

(A:  =  Mfft/2,Mfft/2  +  l,--.,Mfft-l), 


where  /  is  the  frequency  variable. 


4.2  FOCUS  RANGE  VERSUS  MRA  FORMULA 

When  designing  different  sets  of  focus  ranges,  it  is  of  interest  to  keep  the  scalloping  loss 
between  neighboring  range  sets  fixed.  The  width  of  the  aperture  used  for  focusing  determines 
its  range  sensitivity.  Since  the  effective  width  of  the  aperture  seen  by  a  target  located  at 
bearing  0  is  a  function  of  sin^,  and  since  a  curvature  across  an  aperture  (i.e.,  range  sensitiv¬ 
ity)  varies  as  a  function  of  the  square  of  the  aperture  width,  a  focusing  range  that  varies  as 
a  function  of  sin^(MRA),  and  hence  k,  must  be  devised. 

Furthermore,  provided  that  the  quantity  R  —  Xi  cos  6  in  equation  (2)  is  not  significantly 
different  from  R  (i.e.,  R  ^  XiCos6),  which  is  true  at  not-so-oblique  bearings,  then  the  time 
delay  can  be  viewed  approximately  as  a  linear  function  of  inverse  range  -  a  feict  that  must 
also  be  used  in  designing  focusing  range  locations.  Thus,  the  scalloping  loss  can  be  kept  con¬ 
stant  between  range  sets  by  maintaining  a  constant  distance  between  consecutive  focusing 
ranges  in  inverse  range  (i.e.,  =  constant).  This  distance,  however,  while  designed 

to  remain  constant  across  ranges,  must  also  be  made  to  change  as  a  function  of  sin^(MRA), 
for  the  reasons  explained  above. 

Utilizing  this  information,  a  particular  formula  was  designed  as  follows: 


FR=  { 


sin^(MRA) 
-ftmax-TvpR-r 
D  sin2(45) 

u  sin^(135) 
^max  TVpR-r 


45°  <  MRA  <  135° 
MRA  <  45° 
MRA  >  135° 


>  r  =  1  (closest), 2,-  •  • ,  ApR  (farthest), 


(17) 


where  iN/pR  is  the  number  of  focusing  range  sets  desired  and  Rmax  is  the  maximum  value  of 
focusing  range  allowed  for  all  MRA  and  r,  where  r  is  the  range  set  index.  Equation  (17) 
is  illustrated  in  figure  2  for  Rmax  =  18,000  yards  and  iVpR  =  5.  An  infinite  number  of 
realizations  is  possible.  If  the  spacing  between  consecutive  focusing  ranges  is  kept  constant, 
there  remains  one  degree  of  freedom,  which  is  the  location  of  one  of  the  focusing  ranges. 
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Figure  2*  Six  Sets  of  Focusing  Ranges  as  a  Function  of 
Beamformer  MRA,  Expressed  in  Range  and  In¬ 
verse  Range 


4.3  USING  K-uj  VALUES  OUTSIDE  THE  ACOUSTIC  CONE 

To  perform  the  focusing  and/or  interpolation,  points  outside  the  acoustic  cone  (waves 
flowing  below  the  speed  of  sound)  must  be  used.  Such  points  will  produce  a  cos(0)  in  the 
k  ^  6  equation  that  is  greater  than  one  in  absolute  value.  This  can  be  dealt  with,  since  in  the 
equations  needed  to  obtain  the  coefficients,  such  as  the  time-delay  and  the  range- vs- bearing 
equations,  only  the  sin^(^)  is  needed,  which  can  be  replaced  by  1  —  cos^(0). 

4.4  K-u  APPLICATION 

The  K-oj  field  obtained  by  the  two-dimensional  FFT  method  briefly  described  above 
is  the  infinite-range  K-uj,  For  a  given  frequency  (bin  io)  and  look  direction  9o  (bin  ko  = 
Mfftdcos  0o/A),  it  is  of  interest  to  use  different  neighboring  points  (beams)  along  the  k-axis 
(i.e.,  differently  steered  beams  focused  at  infinity)  to  produce  a  differently  focused  beam  at 
the  same  point  (frequency  and  steering  bearing)  of  interest.  So,  for  K-u  point  (io^ko),  it  is 
of  interest  to  use  the  L  infinitely  focused  points  (zq,  ko),  {io,  ko  ±  1),  (zq,  ko  ±  2),...,  selected 
around  (zq,  ko),  but  not  necessarily  symmetric  in  k  about  ko,  in  an  interpolation  scheme  like 
the  one  described  above,  to  replace  the  original  infinitely  focused  point  at  (zq,  A:o)  with  a  new 
one  focused  at  a  prescribed  range. 

The  interpolation  scheme  derived  in  section  3  was  implemented  for  different  values  of  L 
for  a  long  line  array  with  equally  spaced  elements.  Target  frequencies  of  Fs/A  and  Es/16, 
where  Fs  =  3.22  of  the  array  design  frequency  and  MRA  location  ^  61®,  were  selected  so 
that  the  target  is  located  exactly  on  a  K-lo  bin.  The  resultant  beampatterns  —  exact  versus 
interpolated  -  are  overlaid,  showing  the  accuracy  of  the  scheme  as  a  function  of  the  number 
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of  interpolation  points.  The  target  was  swept  in  bearing  at  the  MRA-varying  range  defined 
by  equation  (17)  for  the  two  extreme  range-focusing  cases  of  r  =  6  (farthest)  and  r  =  1 
(closest),  with  A^fr  =  6,  as  shown  in  figure  2.  The  exact-beampattern’s  focusing  range  was 
identical  to  target  range  in  each  case.  The  technique  is  illustrated  in  figure  3.  Simulation 
results  for  different  numbers  of  interpolation  coeflScients  and  frequencies  are  shown  in  figures 
4-11. 


Narrowband  Boampattems  (f=Deslgn.  MRA=6bFR=6750Yds) 


Figure  3.  Illustration  of  MSE  Minimization  Method 
for  Interpolating  27  Beams  Equi-Spaced  in 
Wavenumber  Space  and  Infinitely  Focused  to  a 
Range  Matching  the  Target’s  at  6750  Yards 
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(OP)  WMOd 


Figure  Time-Delay  Focused  (V  =  6, 
^FR  =  6^  Narrowband  Beam- 
pattern  (f  =  Fs/4)  vs  Interpo¬ 
lation  of  Three  Infinitely  Fo¬ 
cused  Ones 


Figure  5.  Time-Delay  Focused  (r  =  6, 
^FR  =  6^  Narrowband  Beam- 
pattern  (f  =  Fs/4)  vs  Interpo¬ 
lation  of  Seven  Infinitely  Fo¬ 
cused  Ones 


Actual  wrsus  interpolated  Narrowbanrl  (F  /16  Hz)  Seampattem  (MaA-6i.47S2'.  r-6} 


Figure  6.  Time-Delay  Focused  (r  =  6, 
^FR  =  6^  Narrowband  Beam- 
pattern  (f  =  Fs/16)  vs  Inter¬ 
polation  of  Three  Infinitely  Fo¬ 
cused  Ones 


Actual  veisus  Interpolated  Narrowtand  (F^/16  Hz)  Beamptftem  (MRA^1.47S2"  r^) 


Figure  7.  Time-Delay  Focused  (r  =  6, 
^FR  =  6^  Narrowband  Beam- 
pattern  (f  =  Fs/16)  vs  Inter¬ 
polation  of  Seven  Infinitely  Fo¬ 
cused  Ones 
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ActuaJ  versus  Interpdaied  Narrowband  {F  /4  Hi)  Beampaflem  {MRA*61.0853",  r»1) 


Figure  8,  Time-Delay  Focused  (V  =  1, 
^FR  =  6^  Narrowband  Beam- 
pattern  (f  =  Fs/4)  vs  Interpo¬ 
lation  of  19  Infinitely  Focused 
Ones 


Actual  versus  tntecpolaiad  Narrowband  {F  ri€  Hz)  Beampanem  (MRA>€t.4752'.  ral) 


Searing  (Oeg) 

Figure  10.  Time-Delay  Focused  (r  =  1, 
iVpR  =  6^  Narrowband  Beam- 
pattern  (f  =  Fs/16)  vs  Inter¬ 
polation  of  13  Infinitely  Fo¬ 
cused  Ones 


Actual  versus  Interpoialed  Narrowband  (F  M  hiz)  Beampatlem  (MRA»6l  .0853'.  r>1 } 


Figure  9.  Time-Delay  Focused  (r  =  1, 
iVpR  —  6)  Narrowband  Beam- 
pattern  (f  =  Fs/4)  vs  Interpo¬ 
lation  of  39  Infinitely  Focused 
Ones 


Actual  versus  Interpoialed  Narrowband  (F^/16  Hz)  Beampanem  (MRA«61 .4752“,  r»1) 


Bearing  (Deg) 

Figure  11.  Time-Delay  Focused  (r  =  I, 
Nfr  —  ^)  Narrowband  Beam- 
pattern  (f  =  Fs/16)  vs  Inter¬ 
polation  of  19  Infinitely  Fo¬ 
cused  Ones 


The  error  between  the  interpolated  infinitely  focused  beampatterns  and  the  actual  beam- 
pattern,  given  by  equation  (7),  is  computed  using  the  optimal  values  of  a  and  /?  for  different 
numbers  of  coefficients  L  and  different  desired  focusing  ranges  i?o-  The  result  is  plotted  for 


15 


two  different  target  frequencies  (F,/4  and  F,/16)  in  figures  12  and  13,  respectively. 


Figure  12.  Error  Between  Actual  and  In- 
terpolated  Beampattems  (In¬ 
finitely  Focused)  vs  Number 
of  Interpolation  Coefficients 
for  Six  Focusing  Ranges  at 
Frequency  Fs/A 


Figure  13.  Error  Between  Actual  and  In¬ 
terpolated  Beampattems  (In¬ 
finitely  Focused)  vs  Number 
of  Interpolation  Coefficients 
for  Six  Focusing  Ranges  at 
Frequency  Fs/16 
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5.  CONCLUSIONS 


As  expected,  the  larger  the  number  of  coefficients  used  to  interpolate,  the  better  the 
resultant  beam  approximation.  In  addition,  it  is  clear  from  figures  4-13  that  a  lower  number 
of  interpolation  coefiicients  is  required  at  lower  frequencies  (under  similar  circumstances) 
and  vice  versa.  Furthermore,  as  shown  in  figures  12  and  13,  the  number  of  interpolation 
coefficients  also  increased  as  the  focusing  range  became  smaller.  These  observations  prompt 
the  need  to  use  a  different  number  of  coefficients  for  different  frequency  bands.  Furthermore, 
rather  than  using  infinitely  focused  beams  to  interpolate  to  the  other  A^fr  ranges,  it  would 
be  more  advantageous  to  use  beams  focused  at  range  iVFR/2  (Afr  even)  and  then  interpolate 
up/down  to  the  farther/closer  desired  focusing  ranges.  As  expected,  however,  if  the  focused 
beams  formed  at  range  N-pn/2  were  originally  formed  from  infinite  focus  range  ones,  then 
errors  in  those  interpolations  will  be  compounded  in  further  interpolations.  It  is  important 
to  note  that  the  error  computed  in  equation  (7)  is  not  the  only  criterion  needed  to  determine 
the  number  of  coefficients.  The  beampatterns,  with  their  sidelobe  level  and  mainlobe  width, 
must  also  be  examined. 

It  must  be  noted  also  that  alternatives  to  the  mean  square  error  minimization  method  for 
beam  interpolation  exist.  The  convolution  method  employed  by  Bernecky  [2],  for  example, 
while  not  as  accurate  as  the  MSE  minimization  method,  is  computationally  much  faster. 
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APPENDIX  A 

VECTOR  AND  MATRIX  DEFINITION  FOR 
JN  ON-SHADING-INCLUDED  COEFFICIENTS 


Sometimes,  it  is  of  interest  to  obtain  focusing-steering  coefficients  that  don’t  include  the 
shading  weights.  It  can  be  shown  that  this  can  be  simply  implemented  by  using  a  different 
definition  than  (5)  and  (6)  for  the  vectors  c  and  s  and  the  matrices  C  and  S.  These  are 
stated  here  without  derivation: 


The  formula  used  to  compute  the  coefficients  is  equation  (9)  in  the  main  text. 


A-1  (A-2  blank) 


APPENDIX  B 

MATLAB  CODE  FOR  COEFFICIENT  COMPUTATIONS 


%  MATLAB  CODE  TO  OBTAIN  BEAM/K  INTERPOLATION/FOCUSING  COEFFICIENTS 

y. 

%  AUTHOR:  Sami  Deeb,  CODE  2123 

function  coeff icients=gen_coeff (N,M_fft ,M,f ,d,c,RO.R»cos_thO,cos_th,shad_embded) 


y.  FUNCTION  WHICH  USES  M  AVAILABLE  LINE-ARRAY  SINGLE-FREQUENCY  BEAMS,  STEERED  TO 
y.  GIVEN  DIRECTIONS,  AND  FOCUSED  AT  GIVEN  RANGES,  MEASURED  AT  ONE  FREQUENCY 
y,  VALUE,  TO  PRODUCE  A  NEW  BEAM  STEERED  TO  A  DESIRED  DIRECTION  AND  FOCUSED  AT 
y.  A  DESIRED  RANGE. 

y. 

y.  N  NUMBER  OF  ARRAY  ELEMENTS 

%  M_fft  THE  SMALLEST  RADIX-2  NUMBER  THAT  IS  LARGER  THAN  THE  NUMBER  OF  ELEMENTS 

y.  f  BEAM/BEAMPATTERN  FREQUENCY 

y.  d  ELEMENT  SEPARATION 

y.  c  SPEED  OF  SOUND 

y.  M  NUMBER  OF  COEFFICIENTS  (ODD) 

%  cos.th  IxM  VECTOR  DENOTING  THE  COSINE  OF  ANGLES  WHERE  THE  AVAILABLE 
y,  BEAMS  ARE  STEERED  TO. 

%  cos_tliO  COSINE  OF  ANGLE  OF  DESIRED  BEAM 

%  R  IxM  VECTOR  REPRESENTING  THE  RANGE  WHERE  THE  AVAILABLE  BEAMS  HAVE 

%  BEEN  FOCUSED  TO. 

%  RO  RANGE  OF  DESIRED  BEAM 

y.  shad.embded  DENOTES  WHETHER  THE  GENERATED  COEFFICIENTS  SHOULD  INCLUDE  THE 
y,  SHADING  OPERATION  (shad_embded=’y' )  OR  NOT  (shad_embded=’n’) 

% 

y,  sphr_td(element_coordinate,cos(th),R,c)  A  FUNCTION  THAT  COMPUTES  THE  TIME 
%  DELAY  OF  A  WAVE  FRONT  TRAVELING  AT  SPEED  c,  INCIDENT  ON  THE  ARRAY  AT 

%  AN  ANGLE  WHOSE  COSINE  IS  cos_th,  AND  FROM  A  RANGE  R,  AT  AN  ARRAY 

%  ELEMENT,  DEFINED  BY  ITS  COORDINATES  element.coordinate . 


shade  =  taylor(N,23,6) ; 
shade  =  shade ( : ) ’ /sum ( shade ) ; 
zero_pad  =  zeros (1, (M_fft-N)/2) ; 
shade  =  [zero_pad, shade, zero_pad] ; 
uniform  =  [zero_pad,ones(l,N) ,zero_pad] ; 
shade  _sq  =  shade.  "'2; 


y.  TAYLOR  SHADING  SL=23  dB,  nbar  =  6 
y.  CONVERT  TO  ROW/NORMALIZE  ~  IxN 
y.  LEFT  AND  RIGHT  ZERO-PADDING  ARRAY 
y.  ZERO-PADDED  lxM_fft  TAYLOR  ARRAY 
y,  ZERO-PADDED  lxM_fft  UNIFORM  ARRAY 
%  IxM.fft 


if  (shad.embded  ==  ’n’)  CC_SS  =  shade_sq;  cc_ss  =  shade_sq;  y,CC_SS:  IxM.fft 


B-1 


else 

end 

CC_SS  =  uniform; 

cc_ss  =  shade; 

‘/.cc_ss:  lxM_fft 

elmnt 

=  d*[0:M_fft-l]-(M_fft-l)*d/2; 

‘/.COMPUTING  ELEMENT 

COORDINATES 

V 

=  zeros (M_f f t, M) ; 

tmp_CC 

CC 

=  zeros (M_f ft, M, M) ;  tmp_SS  = 

=  zeros (M,M);  SS  = 

zeros (M_fft,M,M) ; 
zeros (M,M) ; 

col_ones  =  oiies(M,l),* 

for  ‘/.ELEMENT  TIME  DELAY/PHASE  COMPUTATION  LOOP 

VCm,  :)=2*pi*f*(sphr_td(elinnt(m)  ,cos_thO,RO,c)-  . . . 

sphr_td(elmnt(m) ,cos_th,R,c)) ; 

end 

for  m=l:M_fft  ’/, ELEMENT  LOOP 

~  col_ones*V(m, ;) ;  ‘/.temp  :  MxM  MATRIX  OF  IDENTICAL  ROWS 

temp2  =  temp  -  temp’;  •/.temp2  =  V(m,j)  -  V(m,i) 

tmp_CC(m,:,:)  =  CC_SS(m)*cos (temp2) ; 
tmp_SS(m, : , : )  =  CC_SS(m)*sin(temp2) ; 
end 

‘/.  SUMMING  tmp_CC  AND  tmp.SS  OVER  THE  FIRST  DIMENSION  REPRESENTING  ELEMENTS - 

CC  =  reshape (sum(tmp_CC),M,M);  SS  =  reshape (sum (tmp.SS) ,M,M) ; 

cc  =  (cc_ss*cos(V)) ’ ;  ss  =  (cc_ss*sin(V)) ’ ;  %ss  &  cc:  Mxl 

inv.CC  =  CC\eye(size(CC));  .  ‘/.inv.CC  =  inv(CC); 

a=(CC+SS*(CC\SS))\(cc-SS*(CC\ss)) ;  '/,a=  inv(CC+SS*inv_CC*SS)* (cc-SS*inv_CC*ss)  ; 

b=  -CC\(ss+SS*alpha)  ;  •/,b=  -  inv_CC*(ss+SS*alpha) ; 


coefficients  =  a  +  i*b; 


‘/.Mxl 
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